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ABSTRACT 

This  paper  describes  the  development  and  application  of  a 
combined  detailed  three-dimensional  large  eddy  simulation 
(LES)  and  one-dimensional  analysis  tool  to  predict  and  actively 
control  combustion-driven  dynamic  instabilities  in  gas  turbine 
combustors.  The  integration  of  detailed  finite-rate  kinetics  into 
LES  and  use  of  In-situ  Adaptive  Tabulation  (IS AT)  to 
efficiently  calculate  multi-species  finite-rate  kinetics  in  LES 
along  with  the  use  of  global  kinetics  in  the  one-dimensional 
analysis  tool  was  demonstrated.  The  results  showed  that  LES 
can  be  effectively  used  to  simulate  complex  reacting  flows  in 
gas  turbine  combustors  and  to  identify  regions  of  combustion 
instabilities.  The  results  also  showed  that  the  one-dimensional 
combustor  analysis  with  global  kinetics  can  then  be  used  both  to 
capture  the  combustor  unstable  modes  of  the  predicted  regions 
of  instabilities  and  to  actively  control  these  instabilities.  In 
particular,  the  results  demonstrated  that  by  modulating  the 
primary  fuel  injection  rates  and  the  time-lag  between  the  instant 
of  fuel-air  mixture  injection  and  heat  release,  damping  out  the 
instabilities  may  be  achieved. 

NOMENCLATURE 

a*  = Reference  speed  of  sound 
I)  = Diffusion  coefficient 
K0  = Reaction  rate  constant 
k = Bulk  modulus  of  the  fluid 
L = Combustor  length 
m = Fuel-air  mixture  injection  rate 
Pr,  = Turbulent  Prandtl  number, 


p = Non-dimensional  pressure,  pf  p* 
p*  = Reference  pressure 
p = Combustor  pressure  fluctuations 
Qjh  = Heat  transferred  from/to  the  flameholder 

q o = Heat  of  reaction 

R = Non-dimensional  reaction  rate 

Re'  = Reference  Reynolds  number,  p'a*z/p 

Sc,  = Turbulent  Schmidt  number,  v/D 

T = Non-dimensional  temperature,  t/t* 

T*  = Reference  temperature 

7j  = Non-dimensional  temperature  in  the  ih  numerical  cell 

fh 

T{  = Non-dimensional  temperature  at  the  flameholder  location 
Tign  = Non-dimensional  ignition  temperature 

t = Non-dimensional  time,  t/t* 

t = Reference  time,  if  a 

u = Non-dimensional  axial  velocity,  u/a* 

V = Volume  of  the  fuel  feed  system 
x = Non-dimensional  distance,  x/ L 

y,  = Mass  fraction  of  i,h  species 

Greek  Letters 

a = User  specified  constant  in  Eq.  (7) 

et  = Turbulent  to  laminar  viscosity  ratio,  p,  /p 

(j)  = Equivalence  ratio 
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7 = Specific  heat  ratio,  cp/cv 

/u  = Molecular  viscosity 
fu,  = Turbulent  viscosity 
p = Non-dimensional  density,  p / p 
p * = Reference  density 

ft)  = Non-dimensional  fuel  pulsation  frequency 

Superscripts 

* = Reference  quantities 

- = Dimensional  quantities 

Subscripts 

0 = Total  quantities 

1 INTRODUCTION 

The  Air  Force  Integrated  High  Performance  Turbine 
Engine  Technology  (IHPTET)  combustors  must  operate  free  of 
combustion-driven  dynamic  instabilities  that  could  compromise 
the  structural  integrity  of  high  performance  engines.  To  double 
the  thrust  to  weight  ratio,  as  required  for  future  high 
performance  military  engines,  the  combustor  will  be  required  to 
operate  at  much  higher  overall  design  equivalence  ratios,  P3  and 
T3  than  the  existing  military  aircraft  engines  combustors.  In 
addition,  weight  reductions  in  the  combustor,  diffuser  and  fuel 
injectors  can  be  only  achieved  through  innovative  integration 
and  packaging  of  these  components.  Therefore,  combustor, 
diffuser  and  fuel  injectors  must  be  designed  in  a manner  that 
will  lead  to  instability-free  (or  substantially  damped 
instabilities)  operation.  Active  combustion  control  techniques 
may  also  be  implemented  to  damp  instabilities.  While  the  active 
control  technology  attempts  to  introduce  an  out  of  phase 
disturbance  with  the  combustor  pressure  acoustics,  the  passive 
control  technology  requires  a profound  understanding  of  the 
various  driving  mechanisms  such  as  air  and  fuel  flows 
variations,  unsteady  heat  release  and  their  interaction  that  cause 
combustion  instabilities. 

Most  existing  combustors  design  databases  lack  a 
provision  for  predicting  combustion-driven  instabilities  during 
the  predesign  and  design  phases.  Currently,  a need  exists  to 
predict  and  quantify  combustion  instabilities  in  high 
performance  military  combustors.  Axial,  tangential  and  radial 
instability  modes  may  all  develop  in  the  combustors  that  could 
severely  impact  the  engine  performance  and  its  structural 
integrity.  Various  approaches  are  presently  used  to  predict 
combustion  instabilities.  These  approaches  range  from  one- 
dimensional linear  stability-based  [1-3],  to  one,  two  and  three- 


dimensional non-linear-CFD-based  [4-8],  Mohanraj  et.  al.  [5] 
developed  a one-dimensional  combustor  model  using  a heuristic 
mixing  model  along  with  a semi  empirical  open  loop  active 
controller.  Quinn  and  Paxson  [9]  used  a one-dimensional  model 
to  study  thermo-acoustic  instabilities  in  combustion  systems. 
Other  analytical  models  [10-11]  based  upon  the  unsteady  - 
pressure  wave  equation  in  three-dimensions  were  also 
developed  and  calibrated/anchored  to  experimental  results  under 
controlled  conditions  such  that  extrapolation  to  other  conditions 
can  be  performed. 

In  LES  modeling  of  the  momentum  transport  scales  larger 
than  the  grid  size  are  computed  using  a time-  and  space-accurate 
scheme,  while  the  effect  of  the  unresolved  smaller  scales 
(assumed  to  be  mostly  isotropic)  on  the  resolved  motion  is 
modeled  using  an  eddy  viscosity  based  subgrid  model.  This 
approach  is  acceptable  for  momentum  transport  since  all  the 
energy  containing  scales  are  resolved  and  all  the  unresolved 
scales  (that  primarily  provide  for  dissipation  of  the  energy 
transferred  from  the  large  scales)  can  be  modeled  by  using  an 
eddy  dissipation  subgrid  model.  However,  these  arguments 
cannot  be  extended  to  reacting  flows  since  , for  combustion  to 
occur,  fuel  and  oxidizer  species  must  mix  at  the  molecular  level. 
Since,  this  process  is  dominated  by  the  mixing  process  in  the 
small-scales,  ad  hoc  eddy  diffusivity  concepts  cannot  be  used 
except  under  very  specialized  conditions.  To  deal  with  these 
distinctly  different  modeling  requirements,  a new  subgrid 
mixing  and  combustion  model  has  been  developed  that  allows 
for  proper  resolution  of  the  small  scale  scalar  mixing  and 
combustion  effects  within  the  framework  of  a conventional  LES 
approach. 

The  earlier  studies  [12,  13]  have  established  the  ability  of 
the  LES  model  in  premixed  combustion  and  in  fuel-air  mixing. 
To  reduce  the  computational  cost,  the  past  calculations 
employed  flamelet  models  (for  premixed  combustion)  or 
simulated  fuel -air  mixing  without  detailed  chemical  kinetics. 
However,  for  realistic  simulations  for  the  reacting  flow,  detailed 
finite-rate  kinetics  must  be  used,  especially  if  pollutant 
formation  is  to  be  studied.  As  is  well  known,  the  computational 
effort  involved  when  using  detailed  kinetics  is  so  large  as  to 
make  LES  of  even  a simple  configuration  computationally 
infeasible.  Typically,  global  kinetics  are  employed  to  reduce  the 
computational  cost.  However,  such  kinetics  are  not  able  to  deal 
with  ignition  and  extinction  processes  and  are  also  unable  to 
accurately  predict  pollutant  (Nox,  CO  and  UHC)  formation. 

Due  to  the  limited  data  available  on  combustion-driven 
instability  for  gas  turbine  combustors,  no  consensus  as  to  what 
approach  should  be  used  to  predict  instabilities  during  the 
diffuser/fuel  inject ors/combustor  predesign  and  development 
stages  yet  exists.  The  primary  objective  of  this  study  is  to 
demonstrate  the  application  of  a combined  large  eddy 
simulation  and  an  unsteady  one-dimensional  combustor  model 
approach  for  predicting  the  combustor’s  unstable  modes  and  for 
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investigating  the  effectiveness  of  the  fuel  flow  modulation 
methodology  to  actively  control  and  suppress  the  combustor’s 
instabilities.  The  present  study  is  considered  work-in-progress 
towards  using  LES  along  with  a simpler  analysis  tool  to  predict 
instabilities  at  the  pre-design  and  design  stages  of  gas  turbine 
combustors.  The  paper  is  organized  as  follows.  The  LES  model 
formulation  along  with  its  numerical  implementation  is 
presented  in  the  next  section,  followed  by  sections  on  the 
governing  equations  and  combustion  model  for  the  one- 
dimensional combustor  analysis  tool  and  results.  The  final 
section  is  the  conclusion. 


2 LES  MODEL 

The  compressible  LES  equations  are  obtained  by  Favre- 
filtering  Navier-Stokes  equations  using  a top-hat  filter 
(appropriate  for  finite-volume  schemes).  The  filtering  process 
results  in  terms  in  the  resolved  LES  equations  that  require 
modeling.  The  final  LES  equations  are  avoided  here  for  brevity 
since  they  are  described  elsewhere  [14],  In  general,  the  subgrid 
terms  representing  the  subgrid  stress  tensor,  the  subgrid  heat 
flux,  the  subgrid  viscous  work,  the  subgrid  species  mass  flux, 
and  the  subgrid  enthalpy  flux  all  require  modeling. 

The  subgrid  stresses  in  the  LES  momentum  equations  are 
modeled  using  an  eddy  viscosity,  which  in  turn,  is  modeled  in 
terms  of  the  LES  filter  width  A and  the  subgrid  kinetic  energy 
K***.  A transport  model  equation  for  the  subgrid  kinetic  energy 
is  solved  along  with  the  other  LES  equations.  The  effects  of 
subgrid  turbulence  on  flame  structure  (and  propagation)  in 
premixed  combustion  can  be  quantified  in  terms  of  the  subgrid 
kinetic  energy  and  thus,  it  is  advantageous  to  use  this  type  of 
LES  model  for  reacting  flows.  Another  distinct  advantage  of 
this  approach  is  that  this  model  does  not  assume  equilibrium 
between  subgrid  kinetic  energy  production  and  dissipation 
(implicit  in  algebraic  models)  and  thus,  helps  to  attribute  a 
relaxation  time  associated  with  the  nonequilibrium  in  the 
subgrid  scales. 

To  handle  the  distinctly  different  physics  of  scalar  mixing 
and  chemical  reactions  a new  subgrid  model  has  been 
developed.  Details  of  this  model  have  been  reported  elsewhere 
[12,  13]  and  therefore,  only  summarized  here.  This  model  is 
called  Linear  Eddy  Model  (LEM)  [15]. 


2.1  LEM  MODEL  FORMULATION 

The  details  of  the  LEM  model  have  been  given  elsewhere 
[15]  and  therefore,  are  only  briefly  summarized  here.  LEM  is  a 
stochastic  model,  which  treats  reaction-diffusion  and  turbulent 
convection  separately  but  concurrently.  Reaction-diffusion 
processes  evolve  on  a one-dimensional  (ID)  domain  in  which 
all  the  characteristic  length  scales  in  the  turbulent  field  (from  L 


to  tj)  are  fully  resolved  (6  cells  are  used  to  resolve  ?)).  The 
orientation  of  the  ID  domain  is  in  the  direction  of  the  scalar 
gradient  and  within  this  domain,  the  equations  governing 
constant  pressure,  adiabatic  laminar  flame  propagation  are 
solved.  The  deterministic  simulation  of  the  reaction-diffusion 
processes  can  be  viewed  as  a local  direct  numerical  simulation. 
As  a result,  the  reaction  rate  terms  do  not  require  any  closure. 

Turbulent  stirring  of  the  scalar  field  is  implemented  as 
distinctly  independent  process  that  interrupts  the  deterministic 
evolution  of  the  reaction-diffusion  processes  on  the  ID  domain. 
Stirring  is  implemented  as  stochastic  re-arrangement  events 
called  triplet  maps,  each  of  which  represents  the  action  of  a 
turbulent  eddy  on  the  scalar  fields.  It  has  been  shown  that  this 
mapping  can  capture  correctly  the  physical  increase  in  scalar 
gradient  (without  affecting  the  mean  scalar  concentration)  due 
to  eddy  motion.  Three  parameters  are  needed  to  implement 
these  turbulent  stirring  events:  the  typical  eddy  size,  the  eddy 
location  within  the  ID  domain  and  the  stirring  frequency  (event 
rate).  The  method  to  obtain  these  parameters  is  given  elsewhere 
[12]. 


2.2  NUMERICAL  IMPLEMENTATION  OF  LES  MODEL 

LES  studies  have  been  applied  to  methane-air  flame 
studies  [16-18]  and  to  a generic  combustor  concept.  Methane-air 
flames  similar  to  those  in  the  experiments  [16-18]  are  studied 
here.  In  order  to  obtain  a realistic  chemical  state  over  a wide 
range  of  operating  conditions,  a 15-step,  19-species  skeletal 
mechanism  is  employed.  This  mechanism  (which  included  Nox 
kinetics)  has  been  shown  to  be  quite  accurate  over  a wide  range 
of  equivalence  ratios.  It  is  also  capable  of  predicting  extinction 
and  re-ignition,  which  is  particularly  relevant  here  since  FI 
flame,  is  considered  close  to  extinction  limit.  Table  1 
summarizes  the  chemical  species  and  reactions  in  this  skeletal 
mechanism. 

The  numerical  method  is  the  same  as  in  the  earlier  study 
[12]  and  therefore,  only  briefly  summarized.  To  simulate  a 
stationary  flame,  a moving  observation  window  is  used  that 
translates  with  the  flame  brush  to  maintain  approximately  the 
same  relative  position  between  flame  center  and  observation 
window  (even  though  the  flame  propagates  freely  into  the 
reactants).  All  statistics  are  obtained  relative  to  the  flame  center. 
The  computational  domain  is  chosen  large  enough  to  fully 
capture  this  flame  brush  (typically  6 L).  Earlier  studies  and  the 
present  study  show  that  statistically,  stationary  flames  can  be 
simulated  using  this  approach. 

A fractional  operator  splitting  method  is  used  to  first 
evolve  (via  molecular  diffusion  and  turbulent  stirring)  the  scalar 
equations  (without  the  reaction  terms)  for  a small  time  step  St  to 
obtain  the  thermo-chemical  state  after  mixing.  Then  the  reaction 
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system  is  integrated  over  St  to  get  the  final  scales  and  hence,  the 
computational  cost  for  the  time  integration  can  be  very  high. 
Here,  the  ISAT  approach  is  used  to  efficiently  handle  the  finite 
rate  kinetics. 

In  ISAT,  only  the  accessed  region  of  the  composition 
space,  which  is  a subset  of  the  whole  realizable  region  (the  set 
of  all  possible  combinations  of  compositions  for  a given  number 
of  species),  is  tabulated.  This  tabulation  is  done  as  a part  of  the 
simulation  and  when  the  same  composition  re-occurs,  the  table 
is  searched  and  the  stored  information  is  retrieved  using  fast 
binary  tree  search  algorithm.  Since  only  the  accessible  region  is 
stored,  the  overall  time  required  to  build,  retrieve  and  store 
information  reduces  significantly.  Further  details  of  the  ISAT 
algorithm  is  given  elsewhere  [13]  and  therefore,  avoided  here 
for  brevity. 

To  further  reduce  computational  cost,  the  present  model 
was  parallelized  with  message  passing  (needed  for  the  solution 
of  the  reaction-diffusion  equations)  handled  by  MPI  directives. 
For  the  chemistry  point  problem  no  communication  is  needed. 
Each  processor  builds  a table  for  the  composition  that  occurs 
inside  its  domain  during  the  computation.  This  localizes  the 
ISAT  table  to  each  processor  and  reduces  the  overall  load 
(including  search  and  retrieve  time)  for  each  processor  (as 
opposed  to  building  a single  table  for  all  the  composition  that 
occurs  over  the  whole  of  the  domain).  All  simulations  were 
done  on  a 32  processor  CRAY-T3E  system. 

Computational  efficiency  of  ISAT  is  significant  as 
reported  earlier.  The  size  of  the  table  and  cost  depend  upon  the 
error  tolerance  (which  determined  the  allowable  error  in  each  of 
the  scalar  for  a given  initial  state).  For  higher  accuracy,  this 
parameter  should  be  low,  leading  to  an  increase  in  the  total 
simulations  time  and  storage.  For  a value  of  e,o/  = 0.0008  a 
speed-up  of  around  30  in  the  chemical  update  is  achieved  by 
using  ISAT.  This  is  consistent  with  the  speedup  reported  earlier. 
Table  2 summarizes  the  present  results. 
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The  non-dimensionalization  of  pressure,  p,  temperature,  T, 
density,  p,  and  velocity,  u has  been  obtained  by  using  the 
reference  quantities  p* , T , p* , and  a" . The  coordinate  x has 
been  scaled  by  the  combustor  length  L.  The  time  has  been 
scaled  using  the  sound  wave  transit  time  t*  = hj  a'  . The  heat  of 
combustion  q0  is  assumed  to  be  a constant.  The  ratio  of  specific 
heats  is  denoted  by  y. 

The  non-dimensional  form  of  the  equation  of  state  is  given  as 


p = pT 


(4) 


The  source  vector  in  Eq.  (1)  is  given  as 


3.  GOVERNING  EQUATIONS  FOR  ONE-DIMENSIONAL 
COMBUSTION  INSTABILITY  MODEL 

The  combustor  model  numerically  integrates  the  following 
governing  differential  equations  for  a calorically  perfect  gas. 
The  governing  equations  are  given  below  in  non-dimensional 
form  [9]  as 
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where  the  dependent  vector  w is  given  as 


Equation  (5)  contains  contributions  from  the  reaction,  turbulent 
eddy  diffusion,  wall  viscous  forces,  and  flameholder  heat 
transfer.  The  Reynolds  number,  Re  is  defined  as  p*a  L/p.  The 
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turbulent  viscosity  ratio,  £t  is  defined  as  the  ratio  of  turbulent 

to  molecular  viscosity,  jj/jj.  R is  the  reaction  rate,  defined 
below.  Pr,  and  Sc,  are  the  turbulent  Prandtl  and  Schmidt 
numbers  respectively. 

3.1  COMBUSTION  MODEL  FOR  ONE-DIMENSIONAL 
COMBUSTION  INSTABILITY  MODEL 


order  to  damp  non-physical  oscillations  in  the  vicinity  of  strong 
gradients  such  as  those  brought  about  by  the  combustion  process 
[19-22],  A number  of  grid  independent  studies  were  conducted 
to  ensure  that  physical  instability  would  not  be  damped  [22].  In 
the  current  study,  the  grid  size  was  Ax  = 0.005  mm  for  which  the 
results  were  very  much  grid-independent,  as  will  be  shown  in 
the  results. 


The  combustion  model  employed  in  this  study  was  based 
upon  a single-step  chemistry  reaction  scheme.  The  non- 
dimensional  reaction  rate  has  the  form: 
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where  K(,  is  the  reaction  rate  constant,  Tign  is  the  ignition 
temperature  and  T)  is  the  temperature  in  the  i,h  numerical  cell. 
Due  to  the  one-dimensionality  of  the  combustor  model,  flame 
holding  is  accomplished  by  a 10%  change  of  the  e,  over  the 
combustor  length  from  zero  to  a specified  value. 


3.2  FLAMEHOLDER  HEAT  TRANSFER  FOR  ONE- 
DIMENSIONAL COMBUSTION  INSTABILITY 
MODEL 


3.4  BOUNDARY  CONDITIONS  FOR  ONE- 
DIMENSIONAL COMBUSTION  INSTABILITY 

MODEL 

Boundary  conditions  may  be  imposed  as  either  partially 
opened,  fully  open,  or  choked  inflow  (e.g.  constant  mass  flux) 
ends.  In  either  case  the  model  anticipates  the  flow  direction  and 
applies  appropriate  (e.g.  well  posed)  states.  If  the  flow  is 
outward  from  the  computing  domain,  only  the  static  pressure  is 
imposed.  The  remaining  information,  density,  velocity,  and 
mass  fraction  are  obtained  from  the  interior  of  the  computational 
domain.  If  the  flow  is  inward,  total  pressure  and  temperature, 
and  mass  fraction  are  imposed.  The  remaining  unknown 
quantities  (velocity,  static  temperature  and  pressure)  are 
obtained  through  an  iteration  procedure  by  using  isentropic 
relations  for  static  and  total  temperature  and  pressure.  In  other 
words,  the  velocity  is  first  obtained  by  using  du/d: c = 0 at  the 
exit.  Then,  the  static  temperature  is  obtained  iteratively  through 
the  relation  T0/T  = l + 0.5(/-l)w2 /(yRT).  Finally,  the  static 
pressure  is  obtained  through  p0/pt  relation. 


The  term  Qti,  in  Eq.  (5)  is  represented  by  a simple 
algebraic  expression 

QJh=a(0.9Tign-Tifh)  (7) 

where  a is  a user  specified  constant,  0.9  T,„„  represent  the 
assumed  temperature  of  the  flameholder  to  or  from  which  heat 

may  be  transferred,  and  T/h  is  the  (computed)  gas  temperature 
at  the  flameholder  location.  Note  that  when  this  model  is  used, 
energy  is  not  strictly  conserved  in  the  system.  That  is,  energy 
may  be  brought  into  or  taken  from  the  system  by  this  term.  Tign 
is  assumed  to  be  constant,  since  the  flame  response  time  scale  is 
much  smaller  than  that  of  the  flameholder  wall  temperature. 


3.3  NUMERICAL  METHODOLOGY  FOR  ONE- 
DIMENSIONAL COMBUSTION  INSTABILITY 
MODEL 

The  combustor  model  numerically  integrates  the  above 
equations  of  motion  using  a very  simple,  second-order 
MacCormack  scheme.  Artificial  viscosity  has  been  added  in 


4 RESULTS 

In  this  part  of  the  study,  efforts  were  mainly  focused  on 
demonstrating  the  uncombined  performance  of  the  one- 
dimensional combustor  model  and  the  LES  methodology  in 
suitable  combustor  configurations.  The  objective  here  is  to 
provide  confidence  that  the  two  models  are  capable  of 
predicting  the  reaction  flow  field.  In  sequent  study  Part  II,  the 
information  obtained  from  LES  will  be  used  to  guide  the 
application  of  the  one-dimensional  model  in  a combined 
manner.  Results  obtained  from  the  one-dimensional  combustor 
model  will  be  first  presented,  followed  by  the  results  obtained 
from  LES  respectively. 

4.1  ONE-DIMENSIONAL  RESULTS 

Numerical  results  have  been  obtained  to  demonstrate  the 
applicability  of  the  unsteady  one-dimensional  combustor  model 
to  predict  combustion-driven  instability  in  a generic  gas  turbine 
combustor.  Figure  1 shows  the  schematic  of  the  combustor 
geometry  investigated  in  this  study.  Fully  premixed  fuel-air 
mixture  was  assumed  at  the  inlet  of  the  combustor.  The  fuel 
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type  used  was  JP-8.  Combustor  model  constants  used  in  this 
study  are  given  in  Table  3.  As  indicated  earlier,  two  different 
types  of  exit  boundary  conditions  were  implemented  in  the 
combustor  model,  namely,  full  open  and  choked.  In  this 
preliminary  study,  only  the  open  exit  boundary  condition  was 
used.  The  boundary  condition  imposed  at  the  exit  corresponds  to 
a non-dimensional  pressure  value  of  unity.  In  fact,  the  reference 
pressure  in  this  study  was  assumed  to  be  equal  to  the  back 
pressure.  Due  to  the  one-dimensionality  of  the  combustor 
model,  the  effects  of  the  swirl  velocity  component  was 
accounted  for  by  using  a 1 . 1 % reduction  in  the  static  pressure  at 
the  combustor  inlet. 

The  combustor  model  assumes  that  a constant  fuel-air  flux 
can  be  specified  at  the  inlet  in  the  absence  of  combustor 
pressure  acoustics.  This  condition  is  essentially  a choked 
boundary  in  which  the  effect  of  the  downstream  pressure 
acoustics  are  not  sensed  upstream  of  the  inlet.  However,  a 
choked  inlet  condition  would  lead  to  stable  combustion  as  the 
computations  showed.  The  dependence  of  the  inlet  fuel-air 
mixture  flux  upon  the  combustor  pressure  should,  therefore,  be 
determined  through  a fuel  injection  system  model.  Since  the 
current  study  assumes  the  computational  domain  to  start  at  the 
backward  facing  step,  a lumped-element  parameter  model  for 
the  fuel  injectors/swirlers  was  considered. 


m2  = W] 


pV  dp’ 
k dt 


(8) 


Equation  (8)  expresses  the  dependence  of  the  fuel-air  flow  rate 
at  the  inlet  of  the  combustor,  m2,on  the  combustor  pressure 
oscillations  p' . The  fuel-air  mixture  mass  flow  rate  mx  is  a 
constant.  Equation  (8)  also  implies  that  as  the  combustor 
pressure  oscillation  p' —> 0,  m2  = m , and  a constant  fuel-air 
mixture  flow  rate  is  then  imposed  at  the  combustor  inlet. 

The  location  of  the  fuel  injection  and  the  time  required  for 
the  fuel-air  mixture  to  react  are  very  critical  parameters  for 
either  driving  or  damping  the  instabilities  [6,7],  The  time  delay 
between  the  instant  of  fuel  injection  and  the  instant  of  heat 
release  is  primarily  governed  by  the  mixing  and  chemical 
kinetic  times.  To  account  for  this  time  delay,  various  estimations 
of  this  time  delay  were  made  based  on  the  location  of  the  fuel 
injector  and  chemical  times  for  different  equivalence  ratios 
(i.e.,  0 = 0.5  - 0.8),  initial  temperatures  (i.e.,  Tinle,  = 700  - 900  K) 
and  Damkohler  numbers  (Da  = 1.0  - 100).  The  time  delay  was 
then  used  in  the  combustor  model  to  determine  the  instant  at 
which  the  fuel-air  mixture,  mx , is  injected.  The  resultant 
mixture  injection  was  therefore  shifted  in  time  by  an  amount 


equal  to  the  time  delay  (i.e.,  injection  time  = t±St)  and 
w,  ( t ± St)  is  injected. 

Figures  2 and  3 show  the  combustor  limit  cycle  pressure 
oscillations  and  the  corresponding  heat  release  in  non- 
dimensional  form  and  for  a time  delay  of  1.97xl0"3  s and  at  a 
non-dimensional  location  x = 0.8.  This  value  of  the  time  delay 
was  found  to  produce  the  largest  amplitude  of  the  pressure 
oscillations.  The  cycle  limit  pressure  oscillations  have  a non- 
dimensional  period  of  approximately  0.9  (or  3. 6x1 0‘3  s)  and  a 
corresponding  non-dimensional  frequency  of  1.11  (or  278  Hz). 
It  can  also  be  seen  that  these  strong  oscillations  have  shock 
waves  characteristics,  which  indicate  that  these  are  a strong 
longitudinal  instability  of  the  fundamental  mode  278  Hz.  In 
addition,  a close  examination  of  the  limit  cycle  pressure 
oscillations  and  that  of  the  heat  release  fluctuations  clearly 
exhibits  the  phase  relationship  between  the  heat  release  and  the 
pressure  oscillations.  Pressure  and  heat  release  oscillations  were 
predicted  to  be  in  phase  and  thus  satisfy  Rayleigh’s  criterion  for 
driving  instability. 

Figure  4 shows  the  pressure  field  oscillations  for  a time 
delay  of  3.8xl0‘3  s which  is  equal  to  1.97xl0'3  + half  the  limit 
cycle  period  shown  in  Fig.  2.  It  is  clearly  seen  that  the  strong 
limit  cycle  oscillations  of  Fig.  2 have  been  substantially 
damped,  consistent  with  the  time  lag  theory  of  Crocco  [23],  This 
demonstrates  the  capability  of  the  combustor  model  to  simulate 
non-linear  behavior.  Additional  simulations  were  earned  out  to 
substantiate  the  time  delay  effect  and  the  results  were  all 
consistent.  It  should  be  mentioned  that  when  calculations  were 
made  for  constant  fuel-air  mixture  mass  flow  rate  (time  delay 
and  p'  are  set  equal  to  zero),  a steady-state  combustion  was 
predicted. 

To  further  demonstrate  the  ability  of  the  combustor  model 
to  carry  out  active  control  calculations,  the  combustor  response 
to  fuel  flow  modulation  as  a means  to  actively  control  instability 
was  predicted.  The  fuel  injection  was  pulsed  sinusoidally  at 
different  frequencies  and  amounts.  The  predictions  represent  the 
open  loop  response  of  the  controlled  combustor  as  shown  in 
Fig.l.  Note  here  that  no  phase  shift  and  controller  gain  were 
used  in  the  current  study.  The  inlet  condition  for  the  open  loop 
calculations  was  set  to  be  a constant  fuel-air  mixture  flux  with 
various  time  delays  (i.e.,  non-choked  inlet)  and  the  exit 
condition  was  also  open  with  a specified  non-dimensional 
pressure  of  unity.  The  combustor  length  was  also  kept  the  same 
to  ensure  that  the  combustor  resonant  frequencies  are 
considerably  higher  than  the  fuel  pulsation  frequencies.  The 
amplitude  of  the  amount  of  fuel  pulsed  was  varied  between  2% 
to  10%  of  the  total  primary  fuel  injection  rate  and  the  combustor 
response  was  studied. 

Figures  5-8  show  the  predicted  combustor  response  in 
terms  of  the  non-dimensional  pressure,  temperature,  velocity 
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and  primary  fuel  mass  fraction  for  a time  delay  of  1.97x10'’  s 
and  at  the  non-dimensional  location  x = 0.8.  The  non- 
dimensional  fuel  pulsation  frequency  (O  was  0.6  (i.e.,  <u  = 150 
Hz)  and  the  non-dimensional  magnitude  A was  2%  of  the 
primary  fuel  rate.  The  sinusoidal  fuel  pulsation  was  turned  on  at 
the  start  of  the  computations  (i.e.,  after  one  computational  time 
step  for  convergence  and  code  stability).  The  results  clearly 
show  that  significant  damping  of  the  pressure,  temperature  and 
velocity  oscillations  occurred.  In  particular,  it  can  be  seen  that 
the  pressure  initially  increases  due  to  the  rapid  rise  in  the 
temperature  and  then  it  is  followed  by  a period  over  which  the 
pressure  oscillates  and  finally  the  oscillations  damp  out.  This 
behavior  is  also  reflected  in  temperature  and  velocity  fields. 

Figures  9 and  10  show  the  non-dimensional  pressure  and 
temperature  in  the  combustor  for  A = 10%  of  the  primary  fuel 
rate  and  for  the  same  conditions  as  in  Figs.  5 through  8.  It  can 
be  seen  that  increasing  the  amount  of  pulsed  fuel  does  not 
always  result  in  more  damping.  In  fact,  the  amplitude  of  the 
final  oscillations  appear  to  be  somewhat  greater  than  that  of  the 
previous  case,  for  A = 2%.  However,  the  increase  in  the 
oscillation  magnitude  is  not  significant. 

To  further  study  the  effect  of  fuel  pulsation  on  the 
combustor  response,  the  fuel  pulsation  was  turned  on  after  the 
limit  cycle  pressure  oscillations  were  reached  for  the  case 
presented  in  Figs.  2 and  3.  The  combustor  response  predictions 
are  shown  in  Figs.  11-13.  It  is  clearly  demonstrated  that  if  the 
fuel  is  pulsed  after  the  system  has  reached  a limit  cycle,  the 
decrease  in  the  amplitude  of  the  pressure  oscillations  appears  to 
be  limited  and  never  goes  to  zero  as  the  computations  revealed. 
The  reason  for  this  behavior  may  be  attributed  to  the  inability  of 
pulsed  fuel  to  produce  a secondary  heat  release  oscillation  of 
sufficient  amplitude  and  be  completely  out  of  phase  with  the 
strong  limit  cycle  pressure  oscillations.  This  can  further  be  seen 
in  Fig.  11  where  the  heat  release  fluctuations  are  not  entirely  out 
of  phase  with  the  pressure  oscillations.  It  should  also  be 
indicated  that  a similar  behavior  was  seen  [5],  when  an  active 
control  was  turned  on  after  the  system  has  reached  the  limit 
cycle  behavior.  This  is  indicative  of  the  unsatisfactory 
performance  of  fuel  pulsation  to  control  instability  in  the  limit 
cycle  regime.  Multiple  fuel  pulsation  frequencies  based  on 
multiple  pressure  modes  sensed  in  the  combustor  will  be 
required  to  damp  non-linear  instabilities,  which  will  be 
investigated  in  the  near  future. 

4.2  LES  OF  TRAPPED  VORTEX  COMBUSTOR 

The  TVC  is  a combustor  design  under  investigation  for 
IHPTET  applications.  In  this  concept,  the  combustion  occurs  in 
a vortex  trapped  within  a cavity  and  past  results  suggests  that 
this  can  increase  flame  stability.  Figure  14  shows  schematically, 
the  combustor  simulated  in  the  present  study. 


Various  studies  into  the  application  of  cavity-flow 
interaction  have  been  conducted  in  the  past.  Earlier  [24]  it  was 
noted  that  a cavity-locked  vortex  entrains  very  little  main-flow 
air  and  would  result  in  a low  exchange  of  mass  and  heat 
between  the  cavity  and  the  main  flow,  Since  combustion 
requires  a continuous  supply  of  reactants,  fuel  and  air  must  be 
directly  injected  in  the  cavity  to  sustain  the  burning  processes. 
However,  direct  injection  (mass  addition)  can  disrupt  the  flow 
dynamics,  possibly  resulting  in  cavity  instability.  Previous 
attempts  at  numerical  modeling  of  TVC  [25,  26],  noted  the 
possibility  of  mixing-limited  reactions  in  the  TVC.  For  this 
reason,  fuel/air  mixing  and  its  effect  on  combustion  in  the  TVC 
is  investigated  in  this  study. 

The  model  formulation  and  the  various  closure  models  are 
summarized  in  a recent  paper  [27]  and  therefore,  not  repeated 
here.  In  the  following,  we  briefly  highlight  some  of  the  salient 
features  observed  in  this  study. 

To  investigate  the  fluid  dynamic  mixing  properties  of  the 
TVC,  two  cases  (I  & II)  were  conducted  under  identical 
conditions  except  that  the  air  velocity  at  the  inlet  was  changed 
from  20  m/s  (case  I)  to  40  m/s  (case  II).  Non-reacting  mixing 
was  studied  with  fuel  is  injected  into  the  cavity  with  an 
equivalence  ratio  far  above  stoichiometric  (0  = 4.4).  Thus, 
entrainment  and  mixing  with  the  ambient  air  is  critical  in  order 
to  achieve  efficient  combustion. 

To  increase  mixing  rates,  higher  turbulence  levels  are  also 
desired.  Regions  of  high  mean  and,  perhaps  more  importantly, 
high  RMS  velocities  will  tend  to  enhance  mixing  rates. 
Therefore,  we  investigated  the  effect  of  increasing  in  the 
primary  airflow  rate. 

For  visualization,  the  time  averaged  and  instantaneous 
velocity  vectors  for  cases  I & II  are  shown  in  Figs.  15(a)  & 
15(b),  respectively.  The  higher  annular  flow  rate  (Fig.  15(b)) 
enhances  the  trapped-vortex  strength,  allowing  for  increased 
annular/cavity  flow  interaction.  In  the  time  averaged  results,  it  is 
seen  that  two  vortices  are  present  for  case  I.  In  addition  to  a 
vortex  rotating  in  the  direction  of  the  annular  flow,  the  second 
vortex  rotates  in  the  reverse  direction.  However,  these  vortices 
are  not  observed  in  case  II.  Rather,  only  a single,  large  vortex 
rotates  with  the  flow.  The  instantaneous  images  on  the  other 
hand  do  show  transient  vortices  in  the  shear  layer  and  in  the 
cavity  indicating  a process  of  unsteady  entrainment  of  the 
primary  air  from  the  inlet.  As  discussed  elsewhere  [27]  and 
noted  below,  this  turbulent  mixing  appears  to  be  beneficial  for 
fuel-air  mixing. 

Two  reacting-flow  cases  (IV  & V)  were  also  used  to 
investigate  the  impact  of  fuel/air  mixing  rates  under  reacting 
conditions.  The  primary  equivalence  ratio  of  the  cavity  jets  in 
both  cases  was  again  4.4  and  the  injection  temperature  was  a 
slightly  elevated  value  of  500  K.  As  observed  by  Sturgess  et  al. 
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[28],  the  cavity  flow  entrains  relatively  little  annulus  air  and  will 
tend  to  be  fuel-rich  under  reacting  flow  conditions.  This  was 
confirmed  by  looking  at  the  mean  fuel  and  oxygen  mass 
fractions  [27]  (not  shown  here).  For  the  low  speed  case,  near  the 
injection  ports,  the  RMS  values  are  much  higher  due  to 
fluctuations  in  velocity  and  species  composition.  Far  from  the 
injectors,  however  (near  the  fore-body  wall),  the  oxygen  is 
almost  entirely  consumed  leaving  a fuel-rich  cavity  region.  The 
remaining  fuel  is  finally  consumed  outside  the  cavity  in  either 
the  dump  shear  layer  or  downstream  of  the  after-body, 
depending  on  the  inlet  flow  velocity.  For  the  higher  flow 
velocity,  on  the  other  hand,  more  fuel  is  consumed  in  the  cavity. 
This  indicates  that  the  primary  air  stream  is  entrained  more  with 
increase  in  flow  velocity  and  the  increase  turbulence  enhances 
fine-scale  mixing. 

The  rate  at  which  the  fuel  is  mixed  can  be  viewed  by 
examining  the  species  mixture  fraction,  Z,  as  a function  of 
inflow  velocity.  Shown  in  Figs.  16(a)  and  16(b)  are  the 
instantaneous  and  time  averaged  stoichiometric  mixture  fraction 
surfaces,  Z = Zst,  for  both  reacting  flow  cases  (The 
stoichiometric  surface,  Zst,  has  been  highlighted  by  a single 
black  line,  where  Zst  for  methane  air  is  approximately  0.055).  As 
can  be  seen  in  Fig.  16(b)  (U0  = 40  m/s),  Zst  is  almost  entirely 
contained  inside  the  cavity  region  for  both  the  instantaneous  and 
the  time  averaged  results.  For  the  lower  inflow  velocity,  the 
surface  extends  far  downstream  of  the  cavity  zone. 

Also  shown  in  Figs.  16(a)  and  16(b)  are  the  instantaneous 
and  time-averaged  temperature  color  contours  for  cases  IV  & V. 
Following  the  same  trends  as  Z,  the  peak  temperature  regions 
are  seen  to  shift  from  outside  the  cavity  region  for  a lower 
annular  velocity  to  inside  the  cavity  in  case  V.  In  both  mean  and 
instantaneous  views,  the  cavity  temperature  increases  with 
increases  flow  velocity. 

Comparisons  with  experimental  data  is  possible  for  case 
V.  The  mean  and  RMS  temperature  profiles  reported  by  Hsu  et 
al.  [24]  (obtained  with  CARS  technique)  are  super-imposed  in 
Fig.  17.  The  combustor  conditions  were  nearly  identical  to  those 
simulated  in  case  V,  U0  = 42  m/s,  <|)  = 4.4,  H/d  = 0.59:  however, 
propane  instead  of  methane  was  used  as  the  primary  fuel. 
Despite  this  difference,  the  mean  temperature  trends  should  still 
be  comparable.  A maximum  instantaneous  temperature  of  2025 
K was  obtained  inside  the  cavity  in  case  II  which  is  slightly 
lower  then  the  2150  K reported  by  Hsu  et  al  [24], 

For  both  inflow  velocities,  the  mean  temperature  near  the 
injection  ports  has  not  been  greatly  effected.  However,  away 
from  the  injectors,  the  higher  inflow  velocity  resulted  in  higher 
cavity  temperatures  for  both  mean  and  instantaneous  profiles. 
The  peak  temperature  region  for  the  case  IV  is  not  inside  the 
cavity  but  downstream  in  the  after-body  wake  resulting  in 
longer  thermal  residence  times.  The  longer  residence  times  may 
increase  thermal  Nox  production.  As  would  be  expected  with  the 


higher  inflow  velocity,  the  downstream  temperature  is  lower 
despite  higher  cavity  temperatures.  This  is  mainly  due  to  the 
higher  volume  flow  rate  and  lower  overall  equivalence  ratio. 
The  combination  of  a higher  fuel  consumption  rate  and  higher 
velocity  results  in  a lower  thermal  residence  time.  This  is 
beneficial  to  lower  combustion  pollutant  formation. 

5 CONCLUSION 

Large  eddy  simulation  and  one-dimensional  analysis  were 
carried  out  to  predict  combustion-driven  dynamic  instabilities. 
The  one-dimensional  combustor  model  was  first  used  to  study 
active  control  of  combustion-driven  instabilities  in  a generic  gas 
turbine  combustor.  The  numerical  predictions  clearly 
demonstrated  the  effect  of  the  time  delay  between  the  instant  of 
fuel  injection  and  heat  release.  By  varying  this  time  delay, 
damping  or  driving  of  the  instability  may  occur.  The  predictions 
showed  that  there  existed  a single  time  delay  value  for  the 
combustor  geometry  considered  here  and  for  which  strong 
instability  can  set  out  and  is  manifested  by  a limit  cycle 
behavior.  The  analysis  also  showed  that  these  strong  pressure 
oscillations  can  be  damped  out  by  changing  the  delay  time  by 
some  factor  multiplied  by  the  oscillation  period. 

The  one-dimensional  combustor  model  was  also  used  to 
investigate  the  combustor  response  to  fuel  flow  pulsation  to 
control  instability.  The  results  clearly  demonstrated  that  for  fuel 
pulsation  to  be  effective  in  controlling  instability,  the  fuel 
modulation  must  be  turned  on  at  the  start  of  the  engine. 
Moreover,  the  results  indicated  that  multiple  frequency-based 
fuel  pulsation  might  be  required  to  damp  out  strong  shock  wave 
type  instability  oscillations. 

The  use  of  In-situ  Adaptive  Tabulation  (IS AT)  to 
efficiently  calculate  multi-species  finite-rate  kinetics  is 
demonstrated  in  this  study.  It  is  also  demonstrated  that  the  LEM 
model  can  be  used  to  simulate  premixed  combustion  regardless 
of  the  nature  of  the  flame  structure  (i.e.,  flamelet  or  thin- 
reaction-zones)  without  requiring  any  ad  hoc  fixes.  This  study 
when  combined  with  the  earlier  studies  clearly  suggest  that 
subgrid  combustion  modeling  for  LES  is  a viable  option  even 
when  finite-rate  kinetics  is  to  be  simulated. 

Finally,  it  should  be  mentioned  that  the  combustor  model 
developed  here  is  a very  useful  predictive  tool  that  can  be  used 
during  the  predesign  and  design  stages  to  predict  the  dynamic 
instability  characteristics  of  new  combustor  designs  in  a very 
efficient  manner.  However,  detailed  description  of  the  unsteady 
flow  field  must  be  predicted  in  order  to  use  the  one-dimensional 
model  to  suppress  combustion-driven  dynamic  instabilities. 
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Table  I 15  Step,  19  Species  Mechanism 

H + 0.5O2  = OH 

H2  + O.502  = H + OH 

H02  - 0.502  + OH 

0.502  + #202  = 0H  +H0‘1 

0.502  + 0.502  #2  = H +00 

CHZ  + 00  + 02 #4  = 0.502  + 0#4  + 1-502  #2 

0.502  + 20#3  = #2  + CHi  + 00 

0.5O2  + 0#3  = H + CHiO 

0.5O2  + 0#4  = 0#  + CH3 

0.502  + 00  = 002 

0.5O2  + 02  #6  =0#4  + 0#20 

H + OH  = H20 

H +CH4  + NO  + HCN  = 0.502  + 20#3  + #2 

H + 0.502  + 0#4  +HCN  = 20#3  + #0 

0.5O2  + 0#4  + N H3  + HCN  = H20  + 2 0#3  + #2 


Table  2 ISAT  Performance 


Uoi  — 5E-3 

ffoi  — 8E-3 

Dl 

Total  CPU  time 

16.66  h 

9.46  h 

31.11  h 

Chemistry  time 

0.0108  s 

7.3e-3  s 

0.2278  s 

Speed  up  in 

21 

31 

- 

Chemistry 
Speed  up  in 

2 

3 

_ 

simulation 
No.  of  Records 

40500 

16500 

- 

Table  3 


Fig.  1 A Schematic  of  the  Combustor  and  Control  System. 


Non-Dimensional  Heat  Release  I1  Non-Dimensional  Pressure 


2 Limit  Cycle  Pressure  Oscillations  for  a Time  Delay 
St  =1.97x10°  s. 


Non-Dimensional  Time 


Fig.  3 Non-Dimensional  Limit  Cycle  Heat  Release 
Oscillations  Corresponding  to  St  =1.97xl0'3s. 


Non-Dimensional  Pressure  ” Non-Dimensional  Pressure 
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Non-Dimensional  Time 


Effect  of  Time  delay  on  the  Limit  Cycle  Pressure 
Oscillations  for  a Time  Delay  St  = 1.97xl0'3  + 
1.8x10 3 s. 


Non-Dimensional  Time 


Fig.  5 Effect  of  Fuel  Modulation  on  the  Pressure 
Oscillations  for  Concurrent  Primary  and  Pulsed  Fuel 
Injection. 


Non-Dimensional  Velocity  f Non-Dimensional  Temperature 


la  mi 


co  = 0.6 


to  = 0.6 


A = 0.02 


A = 0.02 


Non-Dimensional  Time 

Effect  of  Fuel  Modulation  on  the  Temperature 
Oscillations  for  Concurrent  Primary  and  Pulsed  Fuel 
Injection. 


Non-Dimensional  Time 

Fig.  8 Effect  of  Fuel  Modulation  on  the  Fuel  Concentration 
Oscillations  for  Concurrent  Primary  and  Pulsed  Fuel 
Injection. 


-0.2 


0 5 10  15  20 

Non-Dimensional  Time 


Fig.  7 Effect  of  Fuel  Modulation  on  the  Velocity 
Oscillations  for  Concurrent  Primary  and  Pulsed  Fuel 
Injection. 


0 5 10  15  20 

Non-Dimensional  Time 

Fig.  9 Effect  of  the  Amount  of  Fuel  Modulation  on  the 
Pressure  Oscillations  for  Concurrent  Primary  and 
Pulsed  Fuel  Injection. 


Non-Dimensional  Pressure  <p  Non-Dimensional  Temperature 
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Non-Dimensional  Time 


Effect  of  the  Amount  of  Fuel  Modulation  on  the  Fig.  12  Effect  of  Fuel  Modulation  on  the  Limit  Cycle  Heat 
Temperature  Oscillations  for  Concurrent  Primary  Release  Oscillations  for  Non-Concurrent  Primary 

and  Pulsed  Fuel  Injection.  and  Pulsed  Fuel  Injection. 


Non-Dimensional  Time  Non-Dimensional  Time 


Fig.  11  Effect  of  Fuel  Modulation  on  the  Limit  Cycle  Fig.  13  Effect  of  Fuel  Modulation  on  the  Limit  Cycle 
Pressure  Oscillations  for  Non-Concurrent  Primary  Temperature  Oscillations  for  Non-Concurrent 

and  Pulsed  Fuel  Injection.  Primary  and  Pulsed  Fuel  Injection 
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Fig.  14  Trapped-Vortex  Geometry  Used  in  This  Numerical 
Study.  Total  Length  (x)  = 285  mm. 


b)  U0  — 40  m/s 


Fig.  15  Non-Reacting  Velocity  Vectors  for  U0  = (a)  20  m/s, 
(b)  40  m/s  (cases  I & II).  Upper  Halves  are 
Instantaneous,  Lower  Halves  Are  Time  Averaged. 


b)  U0  = 40  m/s 


Fig.  16  Reacting  Flow  Temperature  Field  and 
Stoichiometric  Surface  for  U0  = (a)  20  in/s,  (b)  40 
m/s.  Upper  Halves  are  Instantaneous,  Lower  Halves 
Are  Time  Averaged.  Temperature  Color  Contour 
Ranges  Are  300  K (Blue)  and  2000  K (Red). 
Stoichiometric  Mixture  Fraction.  Z = Zs,  s 0.055, 
Shown  with  Black  Line. 


Temperature  (K)  Temperature  (K)  Temperature  (K) 
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Fig.  17  Mean  and  RMS  Temperature  Profiles  at  x = (a)  68, 
(b)  55,  (c)  35  mm  for  U0  = 40  m/s  (Case  V). 

< f > (— ),  fRMS  (-  -),  Teip  (.),  T™s  (A). 
Experimental  results  (elp)  from  Hsu  et  al. 
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PAPER -19,  B.  Sekar 

Question  (J.-P.  Hathout,  USA) 

According  to  the  authors,  pulsing  in  open  loop  causes  the  “unstable”  system  to  become 
“stable”.  Can  the  author  clarify  the  physics  behind  this  behavior?  Can  it  be  inferred 
from  the  dynamics  of  the  model  what  type  of  nonlinearities  are  responsible  for  stability? 

Reply 

The  idea  is,  when  the  system  is  unstable,  it  could  be  damped  by  giving  a frequency  and 
modulating  a portion  of  the  fuel  to  different  amplitude  magnitude.  But  totally  this  cannot 
be  made  to  be  stable.  There  will  still  be  oscillations,  but  either  reduced  or  amplified  de- 
pending on  the  frequency  or  amplitude.  Thus,  you  could  change  the  oscillating  pattern 
and  magnitude  by  different  input  of  either  the  time  delay  or  frequency  or  amplitude.  The 
nonlinearities  to  predict  the  limit  cycle  behavior  comes  from  time  dependent  boundary 
conditions,  time  dependent  formulation  of  the  basic  equations,  and  the  coupling  of  un- 
steady pressure  fluctuations  and  the  fuel-air  oscillations. 


